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In this lecture we compare different QCD-like partition functions with 
bosonic quarks and fermionic quarks at nonzero chemical potential. Although 
it is not a surprise that the ground state properties of a fermionic quantum 
system and a bosonic quantum system are completely different, the behavior of 
partition functions with bosonic quarks does not follow our naive expectation. 
Among other surprises, we find that the partition function with one bosonic 
quark only exists at nonzero chemical potential if a conjugate bosonic quark 
and a conjugate fermionic quark are added to the partition function. 
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1. Introduction 

The QCD phase diagram in the chemical potential temperature plane has 
far reaching phenomenological implications ranging from heavy ion colli- 
sions to the interior of neutron stars. Unfortunately, first principle lattice 
simulations are only possible at zero chemical potential, and our knowledge 
of the phase diagram mainly relies on model calculations (see for exam- 
ple 1, 2). Physically, we know that at zero temperature the baryon density 
is zero below a chemical potential equal to mAr/3. Therefore, in the ther- 
modynamic limit, the QCD free energy and its derivatives, such as for 
example the chiral condensate, do not depend on the chemical potential for 
/i < ttin/3. Since the Dirac operator depends on the chemical potential this 
requires miraculous cancellations in the microscopic theory a problem that 
was coined 3 as the The Silver Blaze Problem. This problem becomes partic- 
ularly manifest in terms of the eigenvalues of the Dirac operator which are 
distributed homogeneously in a strip 4 with a width that increases a function 
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Fig. 1. Possible phases of the QCD partition function at nonzero temperature and 
chemical potential. 

of fi. In this lecture we will mainly focus on the zero temperature axis of 
the phase diagram. To better understand the effect of the baryon chemical 
potential in QCD, we will consider four different partition functions listed 
in Table 1 which were discussed in 5 and 6. Although, the first partition 
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Number of Charged 
Goldstone Modes 
for /i < [i c 


Critical Chemical 
Potential 


(det(D + ^70 + m)) 





/"c = \m N 


(|det(D + ^o + m) | 2 ) 


2 




(det(D+/i7o+m)) 


4 




( | det(D+^7o+m)| 2 ) 


na 


/i c = 



Table 1. Summary of properties of low energy QCD at nonzero chemical poten- 
tial and zero temperature. These partition functions will be denoted by Z Nf—1 , 
Z n= i, Z Nf ~ _1 , Z n= _i, in this order. 



function is physically the most relevant, the other partition functions have 
important applications. Because lattice QCD simulations of full QCD at 
nonzero chemical potential are not possible, one sometimes uses the phase 
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quenched approximation where 

<dct 2 (L> + ^70 + m)) -» (|det(D + ^7o + m)\ 2 ), (1) 

which can be interpreted as a partition function of quarks and conjugate 
quarks 7 . Then Goldstone bosons made out of quarks and conjugate anti- 
quarks have nonzero baryon number resulting in a critical chemical poten- 
tial of m v /2 instead of miy/3. 

The bosonic partition function occurs in the formula for the quenched 
spectral density in the microscopic domain of QCD which is given by 15 

p que„ (ZjM) = (2) 

In a future publication 8 we will consider the expectation value of the 
phase of the fermion determinant given by 

/ 2ie ] _ I dct ( D + M7o + m) \ 

{ 1 \dct(- J D + / i 7 o+m*)/- [6> 

This partition function is not among the above list, but based on our in- 
sights from the bosonic partition function, we will be able to predict the its 
phase diagram. 



2. Gauge Invariance and the Phases of QCD at /x 7^ 

The principle that underlies the independence of the free energy on the 
chemical potential is gauge invariance 9 ' 10 . The Dirac operator can be writ- 
ten as 

D + fj,-y + m = e-^ T (D + m)e^ T , (4) 

which implies that the /x-dependence can be transformed into the boundary 
conditions. A /i-independent free energy is possible in a phase that is not 
sensitive to the boundary conditions. This is the case for /1 < /i c when the 
quarks do not loop around the torus in the time direction. 

Although the partition function Z n —_\ is naively gauge invariant it 
turns out that the regulator of the partition function breaks gauge invari- 
ance so that a ^-independent phase cannot exist. The need for regularization 
is best seen by writing the partition function in terms of eigenvalues 11 

^=-1=/ n^ TT i 2 P(I 2U }) *2 TV ( 5 ) 

where C/C m (e) is the complex plane except two small spheres with radius 
e around ±m. Because of the complex conjugated pole the integral diverges 
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as log e. Instead of this regularization we prefer to regularize the partition 
function as 15 (also known as hermitization 12 ) 

*-=-! = (a*" 1 ( n+ 6 + , D + ™ + m )\. (6) 
\ V —D + jU7o + m e J / 

Since the matrix inside the determinant is Hcrmitian, this partition function 
can be written as a convergent bosonic integral. However, e breaks gauge 
invariance, and for e ^ it is not possible to gauge away the chemical 
potential. We thus find that [i c = in this case. 

Let us now consider the partition function with one bosonic flavor. This 
partition function cannot be written as a convergent bosonic integral and 
therefore cannot be interpreted in terms of bosonic quarks only. The correct 
interpretation is to rewrite this partition function as 6 

z N f =-i = I dct*(£> + M7o + m) 

\ dct(D + ^70 + m)det*(D + ^70 + m) 
and regulate the denominator as in (6). However, contrary to the case of a 
pair of conjugate bosonic quarks, this partition function does not diverge for 
e — > 0, and it is possible to gauge away the chemical potential. In this case 
the free energy will be /x-independent in the thermodynamic limit below the 
lightest particle with nonzero baryon number which is a Goldstone boson 
made out of a bosonic quark and a conjugate bosonic anti-quark. 




3. Low Energy Limit of QCD 

The low-energy limit of the partition functions in Table 1 uniquely follows 
from chiral symmetry and gauge invariance. In Table 2 and Fig. 2 we com- 
pare the bosonic and fermionic (see 7, 10, 13, 14, 15) partition functions 
with a pair of conjugate flavors. 

For Nf = — 1 the partition function is given by the ratio in eq. (7). 
In this case the partition function is finite for vanishing regulator and the 
gauge symmetry (4) is not obstructed. Therefore, we have a /x-independent 
phase for [i < fi c . In this phase the chiral condensate is given by 

(^) N ' — 1 = 1 (Y> k —^— + Y> k — Z k —^—) for M < Mc 

V \ z k +m z k +m z k +mj 

= 2E - S = S. (8) 

For (i > [i c the bosonic contribution to the chiral condensate rotates into a 
pion condensate (see Fig. 3) so that for fi ;§> \i c only the fermionic contri- 
bution remains: 



(tptp) = -S for fi > n c . 



(9) 
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(\det(D + n) + m)\ 2 ) 


(|det(D + M7o +m)\- 2 ) 




no regularization 


regularization 


Goldstone 
Manifold a 


u e 17(2) 


Q e Gl(2)/U(2) 


Chiral 
Lagrangian 




£ki„ = ^TrV^QV^Q- 1 


Covariant 
Derivative 


V U = d U + y[U, B] 


V Q = d Q + ti{Q,B} 




— 2 


Mc = 



Table 2. Comparison of the n — 1 and the n = — 1 partition function. 




Fig. 2. Goldstone spectrum for Z n —\ (left) and Z n — _i (right) 



4. Chiral Symmetry Breaking at y£ and the Dirac 
Spectrum 

In the domain where the kinetic term of the chiral Lagrangian factorizes 
from the partition function, i.e. for /i <C 1/7/ and [z <C l/£ 2 , the quenched 
spectral density satisfies the relation (2). This offers the possibility to test 
our results for Z n —-\ by means of lattice QCD simulations. Calculations 
both with staggered fermions and overlap fermions show an impressive 
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Fig. 3. The /^-dependence of the chiral condensate for Nf 



agreement 16 with (2). 

In the thermodynamic limit the eigenvalues are distributed homoge- 
neously inside the strip 

2p 2 F 2 



|Re(z)| < 



(10) 



Therefore, inside this strip, the quenched chiral condensate goes to zero 
linearly. The behavior of the chiral condensate for full QCD is quite dif- 
ferent. In that case the chiral condensate remains nonzero for m — > 0. On 
the other hand, the eigenvalues still spread out in the complex plane. To 
explain 5 this so called "Silver Blaze Problem" we introduce the "spectral 
density" 



(det(D + ^ 7o + m) £ fc S 2 (z - z k )) 



(11) 



(det(D + jU7o + to)) 

Because of the phase of the fermion determinant p full (z,/i) is in general 
complex. Its microscopic limit is known analytically 17 and can be decom- 
posed as 

p full (z,M) = P quon (^M) + P° sc (^M), (12) 

where p osc is complex with oscillations with a period of 0(1 /V) and an 
amplitude that diverges exponentially with the volume. For V — > oo it 
vanishes outside a region with m < |Re(z)| < |/j, 2 F 2 /S — y. The chiral 
condensate follows the same decomposition 

S full = gquen + s osc_ ^ 
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Fig. 4. The chiral condensate in quenched QCD (S quen ) and in full QCD (S fu11 ) as a 
function of the mass. The support of the Dirac spectrum is in between the vertical lines. 

In Fig. 4 we show the behavior of the different contributions in the thermo- 
dynamic limit. This shows that a nonzero chiral condensate for m — > is 
due to the oscillatory contribution to the spectral density 5 . Therefore these 
oscillations solve the Silver Blaze Problem. 




Fig. 5. Average phase of the fermion determinant as a function of fi. 



5. Conclusions 

The behavior of bosonic partition functions at nonzero chemical poten- 
tial is quite different from what could be expected naively. This surprising 
behavior can be understood from gauging the chemical potential into the 
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boundary conditions. In particular, this shows that the partition function 
with a pair of conjugate bosonic quarks has no /i-indepcndcnt phase. The 
free energy of the theory with one bosonic quark, on the other hand is /Li- 
independent for fjb < However, this partition function only exists as 
a partition function of a pair of conjugate bosonic quarks and a fermionic 
quark with the same mass. Finally, an analysis along the lines of this paper 8 
shows that the expectation value of the phase of the fcrmion determinant, 
eq. (3), behaves as in Fig. 5. 
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